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COHERENT  DETECTION  OF  RADAR  TARGETS 
IN  K-DISTRIBUTED ,  CORRELATED  CLUTTER 


INTRODUCTION 

The  coherent  detection  of  radar  targets  in  a  background  of  correlated  non-Gaussian  clutter  is  a 
difficult  problem  that  has  recently  begun  to  motivate  active  research.  As  a  problem  in  detection 
theory,  this  problem  is  well  understood.  The  optimum  detection  procedure  is  given  by  the  Neyman- 
Pearson  likelihood  ratio.  However,  to  formulate  the  likelihood  ratio,  a  representation,  such  as  a 
closed-form  description,  of  the  multivariate  probability  density  function  (pdf)  is  required.  The  model¬ 
ing  of  this  multivariate  pdf  usually  is  accomplished  by  extending  knowledge  about  the  univariate  pdf 
and  the  covariance  matrix  associated  with  the  process.  In  the  absence  of  knowledge  of  the  complete 
structure  of  the  multivariate  pdf,  this  extension  is  necessarily  arbitrary  and  is  usually  motivated  by  the 
researcher’s  desire  for  mathematical  tractability. 

One  such  extension  of  the  problem  of  detecting  radar  targets  in  non-Gaussian  clutter  has  been 
described  by  Cantrell  [1]  and  by  Farina  et  al.  [2],  The  approach  taken  in  these  studies  was  based  on 
the  transformation  noise  model  described  by  Martinez  et  al.  [3].  This  model  maps  a  given  multivari¬ 
ate  pdf,  usually  Gaussian,  into  a  non-Gaussian  multivariate  pdf  in  such  a  way  that  the  resulting  pdf 
has  the  required  marginal  pdfs  and  the  required  covariance  matrix.  Cantrell  [4]  also  applied  this 
same  technique  to  the  radar  problem  by  starting  with  a  bivariate  complex  Cauchy  pdf  instead  of  a 
Gaussian  pdf.  In  each  of  these  cases,  the  result  was  the  specification,  based  on  the  likelihood  ratio, 
of  a  detection  scheme  that  yields  improved  detection  performance  over  the  performance  obtained  by 
processing  the  same  data  with  a  detector  designed  for  Gaussian  noise.  Although  this  type  of  com¬ 
parison  is  not  a  measure  of  the  optimality  of  the  non-Gaussian  detector,  it  does  reveal  the  benefit  to 
be  gained  by  correctly  identifying  the  clutter  process  as  non-Gaussian  rather  than  assuming  that  it  is 
Gaussian. 


Other  methods  for  constructing  a  multivariate  pdf  with  specified  marginals  and  a  specified 
covariance  matrix  are  available,  and  one  of  these  alternate  methods  is  used  in  this  report.  The 
approach  followed  here  is  motivated  by  physical  arguments  about  the  nature  of  electromagentic 
scattering  from  the  surface  of  the  sea  [5]  and  is  applicable  to  the  problem  of  detecting  radar  targets  in 
sea  clutter.  The  clutter  process  is  modeled  as  a  nonhomogeneous  Gaussian  process  where  the  clutter 
power  level  is  a  random  variable.  This  model  results  in  the  so-called  spherically  invariant  random 
process  (SIRP),  which  has  been  described  in  the  literature  [6-11],  The  problem  of  target  detection  in 
a  SIRP  has  been  the  subject  of  little  investigation.  Yao  [8]  examined  the  case  in  which  the  detection 
threshold  is  set  equal  to  1.  Goldman  [12]  examined  the  conditions  under  which  correlation  receivers 
are  optimal.  Picinbono  and  Vezzosi  [13,14]  examined  the  asymptotic  detection  of  a  weak  signal  in 
the  presence  of  a  real  SIRP  and  found  that  the  asymptotic  detector  is  independent  of  the  pdf  of  the 
unknown  power  level.  Spooner  [15]  examined  the  performance  of  the  optimum  detector  for  the  prob¬ 
lem  in  which  the  form  of  the  pdf  of  the  power  level  is  known  and  the  data  are  comprised  of  indepen¬ 
dent  samples.  None  of  these  results  is  directly  applicable  to  this  study,  although  Picinbono  and  Vez¬ 
zosi  did  derive  the  form  of  the  likelihood  ratio  for  the  problem  considered  here.  This  report  exam¬ 
ines  the  problem  of  radar  target  detection  with  a  small  number  of  samples. 
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In  this  report,  the  multivariate  pdf  that  describes  the  non-Gaussian  clutter  process  is  obtained  by 
averaging  a  multivariate  complex  Gaussian  pdf  with  respect  to  a  random  variable  that  multiplies  the 
normalized  correlation  matrix,  i.e.,  the  variance  or  (in  the  vocabulary  of  sea  clutter  researchers)  the 
power  level.  The  likelihood  ratio  is  then  obtained  in  closed  form.  The  resulting  test  depends  on  the 
unknown  signal  amplitude  and  unknown  initial  phase.  An  approximation  to  this  test  that  is  indepen¬ 
dent  of  the  signal  amplitude  and  initial  phase  is  then  formulated,  and  the  detection  performace  of  this 
approximate  test  is  evaluated  for  some  specific  examples.  This  detection  performance  is  compared  to 
the  performance  of  the  optimum  test  in  which  the  signal  amplitude  and  phase  are  known  and  to  the 
performance  of  the  Gaussian  matched-filter  detector  against  the  same  clutter  data. 

MULTIVARIATE  PROBABILITY  DENSITY  FUNCTION 


A  univariate  pdf  is  assumed  that  describes  the  amplitude  behavior  of  the  complex  clutter  process 
as  well  as  the  normalized  correlation  matrix  of  the  clutter  process.  The  extension  of  this  information 
to  the  multivariate  pdf  proceeds  in  two  steps.  First,  since  the  clutter  process  is  modeled  as  a  complex 
process,  some  assumption  about  the  statistics  of  the  phase  of  the  process  must  be  made.  For  radar 
clutter  modeling,  the  phase  of  the  clutter  is  generally  assumed  to  be  uniformly  distributed.  With  this 
assumption  and  the  assumption  that  the  amplitude  and  phase  of  a  given  sample  are  independent,  the 
univariate  complex  pdf  can  be  immediately  obtained.  This  pdf  must  then  be  extended  to  the  mul¬ 
tivariate  domain. 


The  assumption  of  nonhomogeneous  Gaussian  clutter  automatically  satisfies  the  above-stated 
assumptions.  By  nonhomogeneous  Gaussian  clutter,  we  mean  that  on  any  given  detection  decision, 
the  observation  samples  are  drawn  from  a  population  whose  statistics  are  described  by  the  following 
Gaussian  pdf: 


where 

A 


f(x  I  t) 


1  1 

- exp  — — 

(27r)m  |  A I T7”  2 


is  a  positive  definite  normalized  correlation  matrix, 


(1) 


t  indicates  transpose, 


the  bar  indicates  conjugate, 

m  is  the  dimension  of  the  problem  (i.e.,  number  of  pulses), 


x  is  the  received  signal  vector  of  length  m ,  and 

r  is  the  unknown  power  level  that  varies  from  one  oservation  to  the  next  according  to  a  pdf 
f(r). 


The  resulting  multivariate  pdf  that  describes  the  non-Gaussian  clutter  process  is  obtained  as 

/(■*)  =  fQ  fix  |  T)f{T)dT 
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So 


f  (j)dr 
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The  choice  of  / (r)  influences  the  statistical  behavior  of  this  process.  To  ensure  that  this  statisti¬ 
cal  behavior  matches  the  desired  univariate  amplitude  pdf  g(|jc  |),  /(r)  is  chosen  as  the  solution  to 
the  integral  equation 


Si\x  I)  =  10  £(l*  I  |  T)f(j)dr 


f(r)dT, 


(3) 


where  g  ( |  x  \  \  t)  is  a  Rayleigh  pdf  conditioned  on  the  parameter  r.  This  integral  equation  is  in  gen¬ 
eral  very  difficult  to  solve.  However,  if  g(|jc  |)  is  chosen  to  be  the  so-called  K  distribution  [16], 
which  has  been  shown  to  provide  a  good  fit  to  sea  clutter  amplitude  statistics  [17,  18],  the  solution  to 
this  equation  can  be  found.  It  is  given  by  the  well-known  gamma  pdf, 


fir)  = 


(4) 


where  r(-)  is  a  gamma  function,  rj  fixes  the  mean  of  the  distribution,  and  v  controls  the  skewness. 
Figure  1  shows  curves  that  describe  the  behavior  of  the  gamma  density  for  different  values  of  v. 


POWER  LEVEL,  r 

Fig.  1  —  Power  level  probability  density  function 


With  the  choice  of  the  gamma  pdf  for  the  variation  of  the  power  level  from  one  observation  to 
the  next,  the  univariate  amplitude  pdf  is  given  by  the  K-distribution 


g(|*  I) 


(5) 
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The  multivariate  complex  pdf  of  the  clutter  process  is  obtained  by  solving  the  integral  in  Eq.  (2)  [19]. 
This  is  given  by 


/(•*) 


2 

V 

m 

a/  —  x*A 

V 

v—tn 

n  v—m 

-/  2v  F<A-ir 

rooczxriAl 

l  2  J 

\f  x  A  X 

’  J 

(6) 


where  K„(-)  is  a  modified  Bessel  function  of  the  second  kind  of  order  v.  Figure  2  shows  curves  that 
describe  the  univariate  amplitude  pdf.  The  covariance  matrix  of  the  complex  pdf  is  easily  shown  to 
be  equal  to  rjA. 


Fig.  2  —  K- amplitude  probability  density  function 


NEYMAN-PEARSON  TEST 

The  radar  detection  problem  is  defined  as  a  problem  in  binary  hypothesis  testing  [20],  and  the 
optimum  detection  procedure  is  given  by  the  well-known  Neyman-Pearson  test.  The  problem  is  for¬ 
mulated  as  follows: 


Select  a  complex  vector  sample  y  from  one  of  two  possible  sta¬ 
tistical  populations 

H0.  y  =  n 
H{:  y  —  n  +  s 

where  n  is  a  complex  vector  of  length  m  of  noise  samples  and  s  is  a 
complex  vector  of  length  m  of  signal  samples  that  is  independent  of 
n .  The  desired  test  must  decide  in  an  optimal  manner  if  the  observed 
sample  y  was  caused  by  noise  alone  or  by  signal  plus  noise.  A  choice 
in  favor  of  Hx  is  said  to  be  a  detection  (either  true  or  false). 
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The  Neyman-Pearson  test  is  optimal  in  the  sense  that  for  a  fixed  probability  of  false  alarm,  it 
maximizes  the  probability  of  detection.  The  test  is  given  in  the  form  of  a  likelihood  ratio, 


Hi 


f(y\H0) 


(7) 


where  f(y  \  H0)  represents  the  pdf  of  y ,  given  hypothesis  H0,  and  f(y  \  H{)  represents  the  pdf  of  y , 
given  hypothesis  Hl.  For  the  specific  problem  studied  here,  the  likelihood  ratio  is  given  by 


\NP 


1 - i 

^ —  (y  ~  s)‘ A  *(y  -  s ) 

V 

v  —  m 

y1  A-1y 

V 

v  —  m 

v—m 

V —  (y  -  S)‘A  x(y  -  s) 

V 

V 

iV  v—m 

1 

< 
« |v 

^4 

> 

i 

vT 

i _ i 

(8) 


Here  the  signal  may  be  written  as 


s 


ae 


(9) 


where  a  is  the  signal  amplitude,  <f>s  is  the  initial  signal  phase,  and  s  is  a  complex  steering  vector  that 
represents  the  pulse  to  pulse  phase  shift  of  the  signal  caused  by  the  Doppler  shift.  This  steering  vec¬ 
tor  is  assumed  to  be  known;  in  practice,  a  bank  of  filters  is  built  to  determine  the  Doppler  shift  of  the 
target.  However,  the  signal  amplitude  and  initial  phase  are  are  assumed  to  be  unknown.  As  a  result, 
the  likelihood  ratio  as  given  above  is  not  suitable  as  a  detection  procedure.  The  given  test  must  be 

(1)  accompanied  by  an  estimator  of  the  unknown  quantities, 

(2)  approximated  to  remove  the  dependency  on  the  unknowns,  or 

(3)  averaged  over  an  assumed  statistical  distribution  of  the  unknown  quantity. 

In  this  problem,  step  (3)  is  mathematically  intractable,  and  a  combination  of  steps  (1)  and  (2)  is  used. 

ESTIMATION  OF  UNKNOWN  SIGNAL  AMPLITUDE 

The  estimate  of  the  unknown  signal  amplitude  is  based  on  the  maximum  likelihood  (ML)  estima¬ 
tor  of  the  signal  amplitude  under  the  assumption  that  the  signal  phase  is  known.  This  ML  estimate  is 
discussed  in  Appendix  A.  The  estimate  is  given  by 


a  = 


U'A-'yl 

sc A 


(10) 


Figure  3  is  a  histogram  of  the  typical  accuracy  of  this  estimate  against  simulated  data.  Note  that  the 
estimate  is  essentially  an  implementation  of  the  matched  filter.  It  is  of  interest  to  examine  the 
behavior  of  this  estimate  under  the  hypothesis  of  noise  only.  Appendix  B  shows  that  the  probability 
of  the  estimate  being  greater  than  or  equal  to  a  given  value  (provided  that  the  hypothesis  of  noise 
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Fig.  3  —  Amplitude  estimates;  v  =  1,  17  =  1,  p  =  0.98,  two  pulses 


alone  is  true)  may  be  solved  for  in  closed  form.  As  a  result,  we  may  establish  a  threshold,  say  Ta, 
such  that  the  probability  of  the  amplitude  estimate  is  greater  than  or  equal  to  this  threshold  when  the 
hypothesis  of  noise  alone  is  true  and  is  equal  to  a  desired  probability  of  false  alarm.  If  the  amplitude 
estimate  exceeds  this  threshold,  we  can  declare  a  detection  without  further  processing. 

APPROXIMATION  OF  NEYMAN-PEARSON  TEST  TO  REMOVE  DEPENDENCY 
ON  THE  SIGNAL  PHASE 

To  remove  the  dependency  of  the  likelihood  ratio  on  the  initial  signal  phase,  the  quadratic  form 


(y  -  5)'A  *(y  -  s) 

is  replaced  by 

|yfA_1y  +  s' A~ls  -  2a  I^A-1?  |  |  . 

If  the  quadratic  form  given  above  is  expanded,  the  crux  of  this  approximation  is  seen  to  be  the 
approximation  of  Re(V  A^y)  by  |s'A_1y  | .  This  is  the  implementation  of  the  matched  filter  when 
the  signal  phase  is  uniformly  distributed. 

EVALUATION  OF  DETECTION  PERFORMANCE 

At  this  point,  the  detection  procedure  may  be  evaluated.  To  facilitate  the  numerical  work,  the 
log-likelihood  ratio  is  evaluated  in  each  case.  The  log-likelihood  ratio  yields  exactly  the  same  detec¬ 
tion  performance  as  the  likelihood  ratio.  It  does  however  require  different  detection  threshold  set¬ 
tings.  Figure  4  is  schematic  of  the  detection  procedure.  As  shown,  the  estimate  of  the  signal  ampli¬ 
tude  is  first  formulated.  If  this  estimate  exceeds  a  threshold,  a  detection  is  declared  and  processing 
stops.  If  the  estimate  does  not  exceed  the  threshold,  the  approximate  log-likelihood  ratio  test  is  for¬ 
mulated  and  a  detection  decision  based  on  it  is  made. 
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b) 


DECLARE  TARGET 


DECLARE  NO  TARGET 


I S1  A  1  y  I 

II*  A-1  §  I 


Fig.  4  —  Detection  procedure  for  approximate  detector 


For  the  examples  evaluated,  Pfa  =  10  7,  m  =  two  pulses,  p  (the  correlation  between  the  two 
pulses)  =  0.98,  ij  =  1,  p  =  0.1,  1,5,  and  the  steering  vector  is  given  as 

?  =  (1  +7O  0  +  7 1) , 

which  represents  a  90°  phase  rotation  between  pulses.  Finally,  the  signal-to-clutter  ratio  is  computed 
as 


s/c  —  a2/2r] . 

As  stated  above,  the  threshold  at  the  output  of  the  amplitude  estimator  may  be  obtained  in  closed 
form  (Appendix  B).  Figure  5  shows  curves  describing  this  result  for  the  three  examples  computed 
here. 


A  closed-form  evaluation  of  the  detection  threshold  at  the  output  of  the  log-likelihood  ratio  test 
did  not  prove  possible,  but  the  thresholds  were  determined  by  using  importance-sampling  techniques 
[21,22].  Figure  6  shows  curves  describing  these  results  for  the  three  examples  computed.  To  ensure 
that  the  probability  of  false  alarm  (Pfa)  at  the  output  of  the  approximate  detector  is  essentially  equal 
to  the  desired  value  (1(P7),  the  thresholds  were  chosen  so  that  the  Pfa  caused  by  the  amplitude  esti¬ 
mation  procedure  alone  and  the  Pfa  caused  by  the  evaluation  of  the  log-likelihood  function  alone  were 
each  equal  to  0.5  x  10-7.  This  setting  ensures  the  Pfa  at  the  output  of  the  approximate  detector  is 
<  10”‘7.  Table  1  shows  the  actual  threshold  settings  used  to  compute  the  various  curves. 
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Fig.  5  —  Threshold  settings  for  Ta 
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Fig.  6  —  Detection  thresholds  at  output  of  approximate 
log-likelihood  ratio  test 
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Table  1  —  Threshold  Settings  Used  to  Calculate  Various  Curves 


Approximate  Test,  Xo 


V 

T 

x  a 

T 

0.1 

4.5 

21.0 

1 

1.85 

19.8 

5 

1.16 

18.1 

Gaussian  Matched 
Filter,  Xmf 


V 

T 

0.1 

217 

l 

89.4 

5 

57.1 

Optimum  Neyman-Pearson  Test,  \NP 


s/c 

T 

6.53 

12.9 

4.95 

12.5 

0.51 

12.0 

v  =  0.1 


s/c 

T 

2.10 

12.5 

0.51 

14.0 

-0.73 

14.1 

-2.18 

14.1 

-3.01 

14.0 

s/c 

T 

-9.03 

9.1 

-4.95 

13.5 

-3.01 

13.8 

-1.42 

12.3 

v  =  5 


=  1 


To  evaluate  the  detection  performance,  the  Monte  Carlo  technique  was  used.  The  data  were 
generated  as  follows: 

(1)  a  random  number  was  generated  from  the  distribution  describing  the  power  level  variation: 

(2)  a  pair  of  complex  Gaussian  samples  with  the  desired  correlation  and  the  power  level  from 
step  (1)  was  generated: 

(3)  a  signal  of  a  given  amplitude,  with  random  initial  phase  and  the  desired  phase  rotation 
from  pulse  to  pulse  was  added  to  the  clutter  data;  for  clutter  only,  the  amplitude  was  set  to 
0. 

This  sequence  results  in  two  observations  and  thus  one  detection  decision.  Figures  7  through  9  show 
performance  curves  based  on  1000  repetitions.  The  performance  curves  of  the  optimum  Neyman- 
Pearson  test  and  the  performance  curves  of  the  matched  filter  are  also  shown.  The  performance 
curves  for  the  optimum  test  assume  that  the  signal  amplitude  and  the  initial  phase  are  known.  The 
Monte  Carlo  technique  was  used  to  compute  the  performance  of  the  test  at  several  points;  the  other 
values  of  the  curves  were  interpolated  based  on  these  points.  As  shown  in  Fig.  9,  the  performance  of 
the  optimum  test  for  low  values  of  signal-to-clutter  ratio  was  extrapolated.  This  extrapolation  was 
necessary  since  the  setting  of  the  detection  thresholds  in  this  regime  proved  to  be  prohibitively 
difficult.  As  a  result,  this  particular  curve  is  not  intended  to  be  an  exact  evaluation  of  the  test’s  per¬ 
formance  but  only  a  heuristic  indication  of  this  performance.  The  performance  curves  for  the 
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SIGNAL-TO-CLUTTER  RATIO  (S/C)  PER  PULSE  (dB) 


Fig.  7  —  Detection  performance,  v  =  5 


SIGNAL-TO-CLUTTER  RATIO  (S/C)  PER  PULSE  (dB) 
Fig.  8  —  Detection  performance,  v  =  1 


Fig.  9  —  Detection  performance,  ^  =0.1 
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matched- filter  case  were  also  obtained  by  Monte  Carlo  simulation.  However,  the  detection  thresholds 
for  this  case  were  obtained  in  closed  form,  as  indicated  in  Appendix  B. 

The  curves  indicate  that  the  optimum  test  can  detect  small  signals  in  spiky  clutter  (the  clutter  is 
said  to  get  spikier  as  v  —  0)  much  better  than  the  matched  filter  can,  although  the  performances 
approach  each  other  as  the  signal  strength  increases.  The  approximate  test  yields  performance  that 
lies  between  these  two  cases.  Also,  performance  of  the  approximate  test  and  of  the  matched  filter  are 
seen  to  approach  each  other  as  r  -  oo.  This  occurs  because  the  K  distribution  approaches  the  Ray¬ 
leigh  as  v  -  oo . 


SUMMARY 

A  procedure  for  modeling  a  complex  non-Gaussian  clutter  process  as  a  nonhomogeneous  com¬ 
plex  Gaussian  process  is  described.  The  multivariate  pdf  associated  with  this  clutter  process  is  con¬ 
structed  so  that  it  has  a  desired  marginal  amplitude  pdf  and  a  desired  covariance  matrix.  The  general 
procedure  requires  the  solution  of  an  integral  equation.  For  this  study,  the  K  distribution  is  used  for 
a  marginal  amplitude  pdf,  because  this  pdf  matches  sea  clutter  statistics  and  the  solution  to  the 
integral  equation  for  this  problem  is  known. 

An  approximation  to  the  optimum  Neyman-Pearson  test  is  implemented  and  evaluated.  This 
approximation  is  independent  of  the  signal  amplitude  and  initial  phase.  The  resulting  detection  per¬ 
formance  of  this  approximate  test  is  less  than  that  of  the  optimal  test  but  is  better  than  the  perfor¬ 
mance  of  the  matched  filter  when  applied  against  the  same  clutter  data. 

Because  of  the  arbitrariness  of  constructing  multivariate  pdfs  from  limited  knowledge  of  the 
statistics  of  the  clutter  process,  the  performance  of  the  detector  cannot  be  truly  evaluated  until  it  is 
driven  by  actual  clutter  data.  However,  since  the  clutter  modeling  was  based  on  a  physical  model  of 
the  scattering  of  electromagnetic  waves  from  the  surface  of  the  sea,  this  detector  is  expected  to 
perform  well  against  sea  clutter  data. 

ACKNOWLEDGMENTS 

The  author  thanks  Dr.  Ben  Cantrell  of  NRL  for  suggesting  this  problem  and  for  patiently  moni¬ 
toring  the  progress  of  this  work  during  the  months  that  it  took  to  do  the  research. 

REFERENCES 

1.  B.  Cantrell,  “Radar  Target  Detection  in  Non-Gaussian,  Correlated  Clutter,”  NRL  Report  9015, 

Nov.  1986. 

2.  A.  Farina  et  al.,  “Theory  of  Radar  Detection  in  Coherent  Weibull  Clutter,”  IEE  Proc.  134, 

Part  F(2),  174-190  (1987). 

3.  A.B.  Martinez  et  al.,  “Locally  Optimum  Detection  in  Multivariate  Non-Gaussian  Noise,”  IEEE 

Trans.  Inf.  Theory  IT-30(6),  815-822  (1984). 

4.  B.  Cantrell,  “Detection  of  Signals  in  Non-Gaussian,  Correlated  Noise  Derived  from  Cauchy 

Processes,”  NRL  Report  9086,  Oct.  1987. 

5.  G.R.  Valenzuela  and  M.  Laing,  “On  the  Statistics  of  Sea  Clutter,”  NRL  Report  7349,  Dec. 

1971. 


12 


NRL  REPORT  9130 


6.  A.M.  Vershik,  “Some  Characteristics  Properties  of  Gaussian  Stochastic  Processes,”  Theory  of 
Probability  and  Its  Applications,  9,  353-356  (1964). 

7.  I.F.  Blake  and  J.B.  Thomas,  “On  a  Class  of  Processes  Arising  in  Linear  Estimation  Theory,” 
IEEE  Trans.  Inf.  Theory  IT-14(1),  12-16  (1968). 

8.  K.  Yao,  “A  Representation  Theorem  and  Its  Application  to  Spherically  Invariant  Random 
Processes,”  IEEE  Trans.  Inf.  Theory  IT-19(5),  600-608  (1973). 

9.  B.  Picinbono,  “Spherically  Invariant  and  Compound  Gaussian  Stochastic  Processes,”  IEEE 
Trans.  Inf.  Theory  IT-16(1),  77-79  (1970). 

10.  E.  Conte  and  M.  Longo,  “Characterization  of  Radar  Clutter  as  a  Spherically  Invariant  Random 
Process,”  lEEProc.  134,  Part  F(2),  191-197  (1987). 

11.  H.  Brehm  and'  W.  Stammler,  “Description  and  Generation  of  Spherically  Invariant  Speech 
Model  Signal,”  Signal  Processing  12,  119-141  (1987). 

12.  J.  Goldman,  “Detection  in  the  Presence  of  Spherically  Symmetric  Random  Vectors,”  IEEE 
Trans.  Inf.  Theory  IT-22(1),  52-59  (1976). 

13.  B.  Picinbono  and  G.  Vezzosi,  “Detection  d’un  Signal  Certain  dans  un  Bruit  Non  Stationaire  et 
Non  Gaussien,”  Annales  des  Telecommunications  25,  433-439,  (1970). 

14.  G.  Vezzosi  and  B.  Picinbono,  “Detection  d’un  Signal  Certains  dans  un  Bruit  Spheriquement 
Invariant  et  Caracteristiques  des  Recepteurs,”  Annales  des  Telecommunications  27,  95-110, 
(1972). 

15.  R.L.  Spooner,  “On  the  Detection  of  a  Known  Signal  in  a  non-Gaussian  Noise  Process,”  J. 
Acoust.  Soc.  Am.  44,  141-147,  (1968). 

16.  E.  Jakeman  and  P.N.  Pusey,  “A  Model  for  Non-Ray leighan  Sea  Echo,”  IEEE  Trans.  Antennas 
Propag.  AP-24(6),  806-814  (1976). 

17.  K.D.  Ward,  “Compound  Representation  of  High  Resolution  Sea  Clutter,”  Electron.  Lett. 
17(16),  561-563  (1981). 

18.  K.  D.  Ward,  “A  Radar  Sea  Clutter  Model  and  Its  Application  to  Performance  Assessment,” 
Radar-82,  IEE  Conf.  Pub.  216,  London,  Oct.  1982,  pp.  203-207. 

19.  G.N.  Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions,  2nd  ed.  (Cambridge  University 
Press,  1944),  p.  183. 

20.  M.I.  Skolnik,  Introduction  to  Radar  Systems,  2nd  ed.  (McGraw-Hill  Book  Company,  New 
York,  1980),  p.  376. 

21.  R.L.  Mitchell,  “Importance  Sampling  Applied  to  Simulation  of  False  Alarm  Statistics,”  IEEE 
Trans.  Aerospace  Electron.  Syst.,  AES-17,  15-24  (1981). 

22.  A.  diVito,  G.  G.  Galati,  and  D.  Iovino,  “A  Comparison  of  Variance  Reduction  Techniques  for 
Radar  Simulation,”  Radar-87,  IEE  Conf.  Pub.  281,  London,  Oct  1987,  pp.  510-514. 


13 


Appendix  A 

AMPLITUDE  ESTIMATION 


Under  the  signal  hypothesis,  the  pdf  of  the  observed  vector  y  may  be  written  as 


fly)  = 


i 


exp 


1  (y  -  s)1 A  *(y  -  s) 


(2ir)m  I  A  | 


Tm 


f  < J)d  T , 


where 


(Al) 


s  =  ote^s  . 


The  maximum-likelihood  estimate  of  the  signal  strength  a  occurs  at  ^  =  0.  Assuming  that 

da 

the  partial  differentiation  may  be  taken  inside  the  integral,  we  have 


df(y)  =  1  f°°  f(f)  d 

da  (2x)m  |  A  |  0  r™  da 


1  (y  -  s)‘ A  *(y  -  s) 


exp 


(A2) 


Expanding  the  quadratic  form  and  performing  the  derivative  shows  that  the  above  expression  equals  0 
for  arbitrary  /  (7)  if 


Re(^rA  *y)  —  as '  A  *5  =0. 


(A3) 


Solving  this  equation  yields 


a  —  Re(sA  ly)/s A  *5  . 


(A4) 


for  the  estimate  of  a. 

Implementation  of  this  estimator,  however,  requires  knowledge  of  the  signal  initial  phase.  To 
account  for  our  lack  of  knowledge,  we  implement 

a  =  |  s*  A~'y  |  /s  A ~ls  .  (A5) 
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Appendix  B 

THRESHOLD  SETTING  FOR  THE  MATCHED  FILTER 


In  this  report,  the  form  of  the  matched  filter  considered  is 


xmf  =  I  §  A  y\/a. 


(Bl) 


where 


a  =  1 


for  the  detection  curves 


a  =  s  A  ls  for  the  setting  of  Ta . 


Conditioned  on  the  value  of  r  and  under  the  noise  hypothesis,  the  pdf  of  \mf  is  easily  shown  to  be 

K-  f  " 


/(Xmf)  =  -^exp 


where 


(B2) 


0  = 


A*  1  -  1  A 

s  A  s 


Still  conditioned  on  the  value  of  r,  we  have  for  the  probability  of  a  false  alarm 


Pfa/T  =  exp  <  - 


2/3r 


(B3) 


The  unconditional  probability  of  a  false  alarm  is  then 


pf  =  Jo  exp  1  ~2fr  P™' 


r  "n 
V 


oo 

J. 


r"  1  exp 


v_  _  T__ 
V  T  2/3 r 


K„ 


V0 
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The  thresholds  Ta  were  set  by  letting  (3  =  1/s  A  (i.e.,  a  =  s  A  *s) 

the  output  of  the  matched  filter  detector  were  set  by  letting  /3  =  s  A~ls  (i.e., 
were  checked  and  verified  by  Monte  Carlo  simulation. 


,  and  the  thresholds  at 
a  =1).  These  results 
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